Casimir forces between cylinders and plates 
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We study collective interaction effects that result from the change of free quantum electrodynamic 
field fluctuations by one- and two-dimensional perfect metal structures. The Casimir interactions in 
geometries containing plates and cylinders is explicitly computed using partial wave expansions of 
constrained path integrals. We generalize previously obtained results and provide a more detailed 
description of the technical aspects of the approach We find that the interactions involving 
cylinders have a weak logarithmic dependence on the cylinder radius, reflecting that one-dimensional 
perturbations are marginally relevant in 4D space-time. For geometries containing two cylinders 
and one or two plates, we confirm a previously found non-monotonic dependence of the interaction 
on the object's separations which does not follow from pair-wise summation of two-body forces. 
Qualitatively, this effect is explained in terms of fluctuating charges and currents and their mirror 
images. 



I. INTRODUCTION 

Quantum effects like Casimir forces have become in- 
creasingly important as electronic and mechanical sys- 
tems on the nanometer scale become more prevalent 
0, 9- Now mechanical oscillation modes of quasi one- 
dimensional structures such as nano wires or carbon nan- 
otubes can be probed with high precision Generally, 
the behavior of such systems is influenced by the collec- 
tive nature of fluctuation forces: The total interaction of 
a system of objects or particles cannot be obtained by 
simply adding the forces between all pairs. Instead one 
has to consider also 3-body and higher order interactions 
that become increasingly important with decreasing sep- 
arations between the objects. 

So far Casimir electrodynamic interactions have mostly 
been investigated for two objects: parallel plates @, a 
rectilinear piston @, plate-sphere interaction at asymp- 
totically large distances 0] and for all separations only 
recently Q . A previous letter summarized the results for 
a plate and a parallel cylinder While shape and ge- 
ometry can strongly influence two-body Casimir interac- 
tions, it is also important to understand the consequences 
of the non-additivity of fluctuation forces. In addition, 
the extent to which fluctuations are correlated depends 
on the effective dimensionality of the space that can be 
explored by the fluctuations. Therefore, Casimir interac- 
tions are expected to depend strongly on the codimension 
of the interacting objects [§]. 

In this work we concentrate on two central aspects 
of fluctuation forces: Effects resulting from the non- 
additivity and the particular properties of systems with 



a codimension of the critical value of two. We con- 
sider these problems in the context of interactions be- 
tween cylinders and sidewalls. In previous works we 
have demonstrated that Casimir forces in these geome- 
tries have only a weak logarithmic dependence on the 
cylinder radius [l[ and can be non-monotonic [lol , fll| - 
consequences of codimension and non-addivity. Here we 
employ and extend previously developed methods [l[ to 
obtain the exact interaction between two perfect metal 
cylinders and its modification due to sidewalls from a 
partial wave expansion. These geometries are of recent 
experimental interest since cylinders are easier to hold 
parallel and generate a force that is extensive in its length 

In analogy to the cylinder-plate interaction 
we obtain a weak logarithmic dependence 
on the cylinder radii R a for the interaction 
E ~ -hcL/[d 2 \og(d/Ri)\og(d/R 2 )} between two 
parallel cylinders at asymptotically large distance 
d 3> R- We include higher order partial waves to de- 
scribe the crossover between the asymptotic expression 
and the interaction at very short separations where the 
proximity force approximation (PFA) gives the correct 
zcroth order approximation to the Casimir energy. 
For two cylinders of equal radius R, it has the form 



E 



cyl-cyl 
PFA 
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^hcLy/R/(d- 2R) 5 [Hi- When one or 
two perfectly conducting sidewalls are added to a pair 
of cylinders, the force between the cylinders depends 
non-monotonically on the sidewall separation H [ill ]. 
We compute the interaction between the cylinders and 
the cylinder and the sidewall over a wide range of 
separations by employing the method of images and by 
summing numerically a large number of partial wave 
contributions. The non-monotonic behavior is found to 
result from a competition between force contributions 
from transverse magnetic (TM) and electric (TE) modes 
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which induce opposite image sources. The TE and 
TM forces between two cylinders are monotonically 
increasing and decreasing with the separation of the 
plate, respectively; their sum behaves non-monotonic 
because the slopes arc different. 

The rest of the paper is organized as follows. In the fol- 
lowing Section we describe the methodology of our path 
integral approach and derive the elements of the rele- 
vant matrix operators for cylinders and sidcwalls in a 
partial wave basis. In Section IIIII we obtain analytical 
and numerical results for the forces between cylinders 
and sidewalls. Experimental implications and corrections 
for cylinders of finite length are discussed in Section [TV] 
More technical steps of the calculations are relegated to 
the Appendices. 



II. METHODS 

A. Path integral and partial waves 

We consider geometries that are composed of infinitely 
long cylinders and infinitely extended plates that are ori- 
ented such that the cylinder axes are all parallel and co- 
incide with an in-plane axis of the plates that are also 
parallel to each other. Hence, the geometries have one 
continuous translational symmetry. This allows the elec- 
tromagnetic modes to be split into transversal magnetic 
(TM) modes, described by a scalar field with Dirichlet 
(D) boundary conditions, and transversal electric (TE) 
modes, described by a scalar field obeying Neumann (N) 
boundary conditions. After a Wick rotation to the imag- 
inary frequency (qo) axis, the action for the scalar field 
has the simple form 



for N conditions, where the indices a, f3 label the sur- 



dqo 



|V$| 



(1) 



For the implementation of the boundary conditions and 
the computation of the interaction energy, we employ the 
techniques derived in references [H, [HI, LLJ, [HI] . After in- 
troducing an auxiliary field on each boundary to enforce 
the boundary conditions and integrating out the scalar 
field $, one obtains an effective quadratic action for the 
auxiliary fields with kernel 

Af Q/3 (u, u';q ) = G (s a (u), sp(a'); q Q ) (2) 

for D conditions and 

M a p(u,u';q ) = d naiu) d n[i{ul) Go(s a {u),s (u');q o ) (3) 



faces. Here Go (x, x'; go) 



-? |x-x I/4tt|x - x'| is the 



free space Green's Function and s Q (u) is a vector point- 
ing to the a'th surface parametrized by the coordinate 
vector u which describes a surface in 3D, so that it 
stands for two independent parameters, e.g., for a cylin- 
der u = (xi, <ft) with x\ oriented parallel the cylinder axis 
and </> the azimuthal angle. When wc integrate over the 
auxiliary fields, wc finally obtain the Casimir energy for 
D and N modes at zero temperature (l5l | , 



2tt 



dq Tr logiMM- 1 ). 



(4) 



The total electromagnetic Casimir energy is the sum of 
the energies E D and E N . The force between two objects 
separated by a can be computed by differentiating the 
energy, 



F D/N = -^l°°dqoTrM-^ a M. 



(5) 



The trace runs over the coordinates u and the indices a, 
/?. M^ 1 is the functional inverse of M with all surfaces 
infinitely separated from one another. 

Since the surfaces are static and since the system is 
invariant under translations along the ui-dircction of the 
cylinder axes, it is useful to transform the kernel M to 
momentum space where it is diagonal with respect to 
the momenta qo and q±. If we denote by M(qo,qi) the 
Fourier transformed matrix that is non-diagonal with re- 
spect to the remaining momenta q 2 , q' 2l the energy can 
be expressed in terms of the determinant of that matrix, 



eD/n = hcL J logll M( qo , qi )M-\q , qi ) 



hcL 

47T 



qd q \og\\M{ q )M^{q)\\ 



(6) 



with q = \fqlo + q\ and L the overall extent of the system 
along the Ui-axis. Hence M is a block-diagonal matrix 
of infinite size. The blocks are indexed by a, (3 and the 
momenta q 2 , q' 2 which may take discrete or continuous 
values depending on whether the corresponding surface 
is compact (cylinder) or infinitely extended (plate) along 
the u 2 -direction. The matrix elements are defined by the 
Fourier integrals 



(q2\M a/3 (q)\q' 2 ) = 



G(s^u 2 )-sj(u' 2 );q)e-^^<^^ 



(7) 
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for D modes and 



{q.2\M a p{q)\<h) 



«(«2) 



du2du' 2 
2tt 



(8) 



for N modes where G(x _L ;g) is the xi-Fourier trans- 
formed free Green's function and s^(u 2 ) is the projection 
of s a (u 2 ) onto the x 2 -X3~p\a,ne that is perpendicular to 
the direction of translational invariance. 

After obtaining the matrix elements the determinant 
of MM^ 1 needs to be computed. For n objects M and 
M M,^ 1 obviously have n x n blocks each of which is in- 
dexed by q 2 , q 2 . For two objects we make use of the block 
matrix determinant formula 



Afn M12 
M 2 i M 22 



= \\Mn\\\\M22-M 21 M^M 12 \\ 
= \\M22\W\M11 - M 12 M^M 2l \\ .(9) 



For more than two objects one could recursively reduce 
the size of the matrix to compute its determinant. If the 
objects are infinite flat plates, a simpler approach is to 
employ the method of images which amounts to replace 
the free space Green's function by the Green's function 
for half-spaces or slabs. 



the determinant formula of Eq. §§§ and the fact that the 
off-diagonal matrix elements vanish for d — > 00, we get 



I MM-, 1 1| = ||1 - .\/,,,.l/ 2 , M, ( .l/ 



Si 



(15) 



Hence, the interaction energy of the two cylinders can be 
obtained from Eq. ([6]) and the matrix elements 



{m\M 52 Mz}M 2i M^ 



7 X 



,(1,1) > ]K m+n {qd)ZX(2,2)K n+m ,(qd) 



(16) 



with X — D, N standing for Dirichlct and Neumann 
boundary conditions and the definitions 
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>(i,j) 



Ah i) 



l m (Riq) 
K m ,{R jq ) 



(17) 
(18) 



B. Two cylinders 

We begin by considering two cylinders at center-to- 
center distance d of the general geometry shown in Fig. [I] 
The positions of the surfaces are parametrized in the x 2 - 
X3 -plane by 

s 2 {4>) = (Ri sm Ri cos (f) (10) 
s^ - (</>) = (i?2 sin (J) + d, R 2 cos 0) , (11) 

where the parametrization coordinate here is u 2 = 4>- 
With this parametrization, the matrix elements defined 
in Eqs. ([7]), © can be computed straightforwardly. We 
find (for details see Appendix |A|) 

(m|M 22 |m') = 8 mtml I m (R 1 q)K m (R 1 q) (12) 

(m\M 23 \m') - {-l) m+m 'l m (R 1 q)I m ,(R 2 q)K m _ m ,(qd). 

(13) 

for D boundary conditions. The matrix elements of M 2 % 
for N boundary conditions are obtained through differ- 
entiation of the matrix elements of M23 for D boundary 
conditions with respect to R\ and R 2 . The elements of 
the diagonal blocks for N boundary conditions are given 
by 



(m|M 22 |m') = 5 m , m ,q 2 l' m {Riq)K' m (Riq). 



(14) 



The matrix elements of M33 are given by Eqs. (fT2|) and 
(THI) with i?i replaced by R 2 and M 32 = M 23 . Using 



H 



H 



S2 



d 



FIG. 1: Surfaces and parameters used for computing matrix 
elements. Surface 53 is displaced by d to the right of surface 
S2. Surface S3 reflected at the plate yields surface 53, r. 



C. Cylinders and plates: Method of images 

Now wc add infinite plates to the geometry, see Fig. [1] 
The effective number of surfaces can be kept the same 
when adding the plates to the system if instead of the free 
space Green's function, Go, modified Green's functions 
are used, which obey Dirichlct or Neumann boundary 
conditions at the plates. For one plate the half-space 
Green's function is 

G^p(x,x';<? ) = G (x,x';<jo) - s x G (x, x' fl ; q ) (19) 

where x^ = (xi, x 2 , — x 3 + 2H) is the reflection of x' 
at the plane at X3 = H and s D = +1, s 



N 



■1. For 
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two plates infinitely many images have to be used since 
successive reflections at the plates generate a series of 
images with increasing separation from the plates. For 
this slab geometry, the Green's function can be written 
as the series 

G fp( x , x ';<?o) = G (x,x';q ) 
-s A [G (x,x' fll ;go) + G (x,x^ 2 ;<7o)] 
+ [Go(x,x^ 1H2 ;(7o) + Go(x,x^ 2 j ?1 ;< ?0 )] (20) 
-s A [G (x,x' fll M2, R1 ;qo) +Go(x,x' R2yRliR2 ;qo)} 



where x 



Ra,Rp,... 



is obtained from x' by a sequence of 
reflections at plate a, (3, . . . The theory developed for the 
free space at the beginning of this Section applies also to 
half-space and slab geometry. The interaction between 
a set of cylinders in the presence of one or two parallel 
plates can be obtained again from Eq. ([6]) ; one only needs 
to substitute the half-space or slab Green's functions for 
G in Eqs. (0, ((HJ and the obtain cylinder-cylinder matrix 
elements. When using these image Green's functions, 
their appropriate Fourier transformed matrices shall be 
labeled M lp and M 2p . 

Below, we shall consider the two cylinders S 2 and S3, 
see Fig. [TJ The surfaces of the reflected cylinders have 
the paramctrization 



s 2,r(<P) = ( R i sin</>,±i?i 
s Ir(0) = (R2 sin M 



cos <fi + S) 
d, ±i? 2 cos 



4> + 5), 



(21) 
(22) 



where S is a distance normal to the plates. Upon one 
reflection 6 is 2H and the cylinder surface orientation 
is reversed, hence, the minus sign is chosen in Eqs. [21] 
and[22l For two plates the reflected cylinder is reflected 
again, so that the plus sign must be chosen for the orien- 
tation, and S = —2Hi — 2H 2 (if the first plate is loacted 
a distance Hi above and the second plate a distance H 2 
below the cylinder, see Fig. [5]). 

The corresponding matrix elements of Eqs. ([7]), © for 
D modes are 

{m\M 2{2>R) \m') = (Tl) m ' I m {Riq)I m >(Riq)K mTm ,(q5). 

(23) 



{m\M 2{m \m') = (Tl) 

x /,„ (Riq)I m , (R 2 q)K mTm > (q\/ S 2 + d 2 ). 

The elements of M3(z,r) follow from those of M 2 (2,R) by 
replacing Ri by R 2 and the elements of M3( 2j _r) are given 
by those of M 2 i^ R ^ with R\ and R 2 interchanged and d 

replaced by — d. For N modes, the elements of M 2 (3 : r) 
are obtained by differentiation of the D elements with 
respect to i?i and i? 2 . The elements of M 2 (2,R) for N 
modes are given by 

(m\M 2(2tR) \m') = {Tl) m ' q 2 I' m {Riq)I' m i(Riq)K mT r,Aq5) ■ 

(24) 



III. INTERACTION ENERGIES 



A. Two cylinders 
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d 



FIG. 2: Cylinder-cylinder geometry. 

We consider two cylinders of radii R\ and R 2 with 
center-to-center separation d, see Fig. [2j For this ge- 
ometry the interaction energy is obtained from Eqs. ([6|) 
and (Tl~5|) with the matrix elements of Eq. (Ti"6|) . For 
large separations d ^> R\, R 2 , the asymptotic behavior 
of the energy is determined by the matrix elements for 
m = ml = for D modes and m = m' = 0, ±1 for N 
modes. Taking the determinant of the matrix that con- 
sists only of these matrix elements and integrating over 
q yields straightforwardly the attractive interaction en- 
ergies 



E D = — 



hcL 



1 



d 2 87rlog(d/i?i)log(d/i? 2 ) 
1 1 



E 



N 



1 



-hcL 



log(d/i?i) 

7 R\R\ 
5tt d e ' 



log(d/fl 2 



(25) 



The asymptotic interaction is dominated by the contri- 
bution from TM (D) modes that vanishes for R a — * 
only logarithmically. 

For arbitrary separations higher order partial waves 
have to be considered. The number of partial waves has 
to be increased with decreasing separation. A numeri- 
cal evaluation of the determinant and the g-integration 
can be performed easily and reveals an exponentially fast 
convergence of the energy in the truncation order for the 
partial waves. Down to small surface-to-surface separa- 
tions of (d — 2R)/R = 0.1 we find that m = 40 partial 
waves are sufficient to obtain precise results for the en- 
ergy. The corresponding result for the energies of two 
cylinders of equal radius is shown in Fig. [3] Notice that 
the minimum in the curve for the total electromagnetic 
energy results from the scaling by the PFA estimate of 
the energy. The total energy is monotonic and the force 
attractive at all separations. 



B. Cylinder-plate geometry 

The simplest geometry to which the method of images 
can be applied is composed of an infinite plate and a par- 
allel cylinder of radius R, see Fig. 2J The Casimir energy 
for this geometry has been computed in Ref. [l[ from the 
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(d-2R)/R 
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The total electromagnetic Casimir interaction is again 
dominated by the contribution from the D mode with 
m = which depends only logarithmically on the cylin- 
der radius. The interaction at all separations follows, as 
in the case of two cylinders, from a numerical computa- 
tion of the determinant of Eq. (|26|) and integration over 
q. The result is shown in Fig. [5] 



H 



FIG. 4: Cylinder-plate geometry. 



FIG. 3: Casimir energy for two cylinders of equal radius R as 
a function of surface-to-surface distance d — 2R (normalized 
by the radius). The energy is divided by the PFA estimate 
EpFA yl f°r the energy given in the introduction. The solid 
curves show our numerical results; the dashed lines represent 
the asymptotic results of Eqs. 1251 The 1/log corrections to 
the leading order result for TM modes cause very slow con- 
vergence. 



matrix elements of Eqs. ([7]), ©for the plate and the 
cylinder. Here, we employ the method of images so that 
we have to consider only one surface, the cylinder, that is 
placed into the half-space that is bounded by the plate. 
Hence, we substitute G by G\ p in Eq. ([7]), ([5]) together 
with the parametrizations of Eq. (jlip and of Eq. ([22]) with 
the + sign, R 2 = R and 5 = 2H. The resulting M\ v is 
simply one block for the cylinder since there is effectively 
only one surface, and it equals M33 — s x M^^y When 
the plate is moved to infinite separation from the cylin 



der, there is only the free cylinder, so 
The determinant of Eq. ^ can now be written as 



M33 1 . 



\M lp M 



lp.oo I 



1 - s x M mR) M. 



(26) 



The matrix elements of M^j^M^ are given for D 
modes by 

(m\M 3i3R) M^\m') = M^L Km+m ,{2Hq) (27) 



and for N modes by 
{m\M mR) M^\m') 



K m ,{Rq) 

urn 



(2Hq) . (28) 



This result in combination with Eq. ([6]) is identical to 
the one given in Eqs. (5)-(8) of Ref. The asymptotic 
expression for the attractive interaction energy at H R 
reads 



E D 



hcL 



H 2 l6wlog{H/R)' 
5 R 2 



(29) 
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FIG. 5: Casimir energy for one cylinder of radius R parallel to 
one plate as a function of the surface-to-surface distance H—R 
(normalized by the radius). The energy is divided by the PFA 
estimate E^'^ tc of energy given in the Introduction. The 
solid curves reflect our numerical results; the dashed lines 
represent the asymptotic results of Eqs. 1291 Convergence for 
the TM (Dirichlet) energy to the asymptotic result is very 
slow because of 1/log corrections. 



C. Two cylinders, parallel to plate(s) 

The image technique lends itself to studying multi- 
body interactions involving plates and cylinders. We 
consider the geometry shown in Fig. [S] with two cylinders 
that are placed parallel to one or in-between two parallel 
plates. Rodriguez et al. [ic| studied a similar geome- 
try consisting of two metal squares between two paral- 
lel metal sidewalls by computing numerically the mean 
stress tensor and observed that the force between the two 
squares changes non-monotonically when the two plates 
arc pulled away. In previous work, we applied the stress 
tensor method and the path integral approach presented 
above to the geometry of Fig. [H] and found again a non- 
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monotonic dependence of the force between the cylinders 
on the separation between the plates Here we pro- 
vide the technical details of the path integral approach 
and the method of images employed in the latter work. 




FIG. 6: Two cylinders parallel to one plate or sandwiched 
between two parallel plates. 

To study the interaction between the cylinders and 
plates it turned out to be more convenient to compute 
directly the forces between the objects as defined by 
Eq. ([5]). For the force between the two cylinders we set 
a = dk.2 and for the force between the cylinders and a 
plate we choose a = -H1X3. The matrix M of Eq. (|5|) can 
be constructed by the method of images as follows. We 
consider the two cylinders as the surfaces that are placed 
either inside a half-space or a slab so that the matrix el- 
ements of M are given by Eqs. {7}, ([8]) with G replaced 
by G* p or G$ p of Eqs. flT^), (|2"0"]l. respectively. We shall 

refer to the corresponding 2x2 block matrices M\ v and 
M2p matrices as M np which is of the form 

= ft" 22 f " 23 ) • (30) 

In the two-body cylinder case without sidewalls the M aa 
matrix blocks describing the self-interaction were diago- 
nal in m. Here, the self-interaction blocks M nPjQ , Q con- 
tain image information and are not diagonal. The ma- 
trix elements of M np are constructed from the elements 
of Eqs. (|2"3")l . (pM|) . Their explicit form and a formula for 
the inversion of M np are given in Appendix [Bj The inte- 
grand of Eq. ([5]) can be straightforwardly computed by 
truncating the matrix M at a finite partial wave order 
m and performing the matrix multiplication and trace 
in Fourier space. Including up to m = 35 partial waves, 
we obtain for the Casimir force between two cylinders of 
equal radii in the presence of one or two sidewalls the 
results shown in Fig. [71 In this figure the force at a 
fixed surface-to-surface distance d — 2R = 2R between 
the cylinders is plotted as a function of the relative sep- 
aration (H — R) / R between the plate and cylinder sur- 
faces. Two interesting features can be observed. First, 
the attractive total force varies non-monotonically with 
H: Decreasing for small H and then increasing towards 
the asymptotic limit between two isolated cylinders for 
large H, cf. Eq. (|25p . The extremum for the one-sidewall 



case occurs at H — R ss 0.27i?, and for the two-sidewall 
case is at H — R « 0.46i?. Second, the total force for the 
two-sidewall case in the proximity limit H = R is larger 
than for H/R — * 00. As might be expected, the in- 
dependence for one sidewall is weaker than for two side- 
walls, and the effects of the two sidewalls are not addi- 
tive: not only is the difference from the H — > 00 force not 
doubled for two sidewalls compared to one, but the two 
curves actually intersect at a separation of H/R = 1.13. 




(H-R)/R 



FIG. 7: Electromagnetic Casimir force per unit length be- 
tween two cylinders for the geometry of Fig. [6] with Hi = 
H2 = H and Ri = R2 = R vs. the ratio of sidewall 
separation to cylinder radius (H — R)/R, at fixed distance 
(d — 2R)/R = 2 between the cylinders, normalized by the to- 
tal PFA force per unit length between two isolated cylinders 
[•Fpfa = f (fic7r 3 /1920) ^R/(d - 2R) 7 [O]]. The force is at- 
tractive. The solid lines refer to the case with one sidewall, 
while dashed lines depict the results for two sidewalls. Also 
shown are the individual TE (blue) and TM (red) forces. 

A simple generic argument for the non-monotonic side- 
wall effect has been given in Ref. [ill It arises from a 
competition between the force from TE and TM modes 
as demonstrated by the results in Fig. [7] An intuitive 
perspective for the qualitatively different behavior of the 
TE and TM force as a function of the sidewall distance 
is obtained from the method of images. For the D modes 
(TM polarization) the Green's function of Eq. p9[) is 
obtained by subtracting the contribution from the im- 
age so that the image sources have opposite signs. Any 
configuration of fluctuating TM charges on one cylinder 
is thus screened by images, more so as H is decreased, 
reducing the force on the fluctuating charges of the sec- 
ond cylinder. This is similar to the effect of a nearby 
grounded plate on the force between two opposite elec- 
trostatic charges. Since the reduction in force is present 
for every charge configuration, it is there also for the aver- 
age over all configurations, accounting for the variations 
of the TM curves in Fig. [7J 

By contrast, the N modes (TE polarization) require 
image sources of the same sign as demonstrated by the 
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half-space Green's function of Eq. (jl9|) . The total force 
between fluctuating sources on the cylinders is now larger 
and increases as the plate separation H is reduced. (An 
analogous additive effect occurs for the classical force be- 
tween current loops near a conducting plane.) Note, how- 
ever, that while for each fluctuating source configuration, 
the effect of images is additive, this is not the case for 
the average over all configurations. More precisely, the 
effect of an image source on the Green's function is not 
additive because of feedback effects: the image currents 
change the surface current distribution, which changes 
the image, and so forth. For example, the net effect of 
the plate on the Casimir TE force is not to double the 
force as H — > R. The increase is in fact larger than two 
due to the correlated fluctuations. 

In Fig. [51 we show the total force between the cylinders 
vs. the sidewall separation H/R for a variety of differ- 
ent values of the cylinder separation d/R in the presence 
of a single sidewall. As we vary d/R the depth of the 
minimum in the force changes, see Fig. [5] The sepa- 
ration (d — 2R) /R = 2 from Fig. [7] seems to achieve the 
largest value of non-monotonicity. For larger or smaller d 
the degree of non-monotonicity (defined as the difference 
between the minimum force and the force in the limit 
H — > 0) decreases. For small d, the force approaches to 
PFA estimate. For large d, the TM (Dirichlet) force dom- 
inates except when the cylinders arc sufficiently close to 
the metal plate when it is reduced enough by its image 
cylinder that the TE (Neumann) force takes over. 




(H-R)/R 



FIG. 8: Casimir force per unit length between two cylinders 
of equal radius vs. the ratio of sidewall separation to cylinder 
radius (H — R)/R (for one plate), normalized by the total 
PFA force per unit length between two isolated cylinders for 
various (d — 2R)/R. The non-monotonic effect appears to 
become weaker as (d — 2R)/R is moved away from 2. The 
vertical lines indicate the sidewall separations used in Fig. [9] 

While the above arguments explain the competition 
between TE and TM forces, they do not show that the 
sum of these competing forces is nonmonotonic. For ex- 
ample, if the TE and TM variations with H were equal 
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FIG. 9: Total Casimir force between two cylinders of equal 
radius R and a sidewall vs. the ratio of cylinder surface sepa- 
ration to cylinder radius (d — 2R) / R, normalized by the total 
PFA force per unit length between a cylinder and a plate 
Fppa = §(fic7r 3 /960) y/R/2{d - 2R) 7 for plate separa- 
tions of H — R = 0.67?, R, and 2R. Note that the normaliza- 
tion is different from the cylinder-cylinder PFA in the previous 
figures. 



and opposite, they would cancel with no net dependence 
on H . That this is not the case can be checked by ex- 
amining the two limits H — R 3> R and H — R -C R. 
In order to simplify the analysis, we assume that there 
is one sidewall and that in both limits the two cylinders 
have a large separation d ^ H — R, R. Then in the case 
H — R ^> R the cylinders and their images are separated 
by a distance that is large compared to R so that the 
forces are dominated by the lowest partial waves, s-wave 
for TM and both s- and p-wave for TE modes PJ. The 
former is stronger and dominates the asymptotic force 
for which we obtain 



Fd 
L 



Ahc 



H 4 



7T d 7 \ri 2 (R/Hy 



(31) 



confirming the reduced net force as the cylinders get 
closer to the plate. 

It is instructive to justify the scaling of the force in 
Eq. (|3ip with H and d from simple physical arguments. 
While the logarithmic dependence on R could have been 
anticipated the H 4 scaling is a remarkable conse- 
quence of the multi-body effect. For TM modes, the field 
obeys D boundary conditions so that the mirror source 
has opposite sign. Therefore, each cylinder and its mir- 
ror image can be considered as a dipolc of size ~ H. 
The interaction of the two dipoles should scale as the 
interaction between two cylinders of size ~ H with Neu- 
mann boundary conditions. For d H the force for the 
latter problem scales as ~ H A /d 7 , explaining the above 
result [lj, up to the logarithm. 

In the opposite limit H — R <C R, the image cylinder is 
very close to the original cylinder so that the interaction 
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involves partial waves of high order. In an attempt to 
analytically understand this limit, we performed certain 
asymptotic calculations reported in Appendix [Cl Since 
the relevance of the conclusions to the regime studied 
numerically (Fig. [7]) is not clear, we shall not further ex- 
plore the implications of these results. We note, however, 
that the numerical results in the regime of Fig. [7]support 
the dominance and fast decay of TE modes for this range 
of parameters, as justfification for the non-monotonic be- 
havior. 

Thus far we have discussed the variation of the force 
between the two cylinders with the sidewall separation. 
We found that the force is not monotonic in H . This 
also implies that the force between the cylinders and 
the sidewalls is not monotonic in d. A nonmonotonic 
force F X2 between the cylinders means that there is a 
value of H where dF X2 /dH = 0. Since the force is 
the derivative of the energy, F X2 = —d£/dd, at this 
point d 2 £ /dddH = 0. These two derivatives, of course, 
can be interchanged to yield d(d£ / dH) / dd = 0. But 
this means that dF X3 /dd = at the same point, where 
F X3 = —d£/dH is the force between the cylinders and 
the sidewall. This cylinders-sidewall force is plotted in 
Fig. [9] as a function of d/R for various values of H/R and 
clearly is non-monotonic in d/R. The non-monotonicity 
is smaller which is not surprising since the effect of a 
small cylinder on the force between two bodies is smaller 
than the effect of an infinite plate. 



IV. DISCUSSION 

In previous research, unusual Casimir force phenom- 
ena were sought by considering parallel plates with ex- 
otic materials: for example, repulsive forces were pre- 
dicted using magnetic conductors (l7j . combinations of 
different dielectrics [l8|, fluids between the plates 
and even negative-index media with gain [20|. A differ- 
ent approach is to use idealized materials such as per- 
fect conductors with more complicated geometries: as 
illustrated in this and previous [1 01 ] work, surprising non- 
monotonic (attractive) effects can arise by considering as 
few as three objects. These effects arise from the collec- 
tive properties of fluctuation forces and cannot emerge in 
a system of particles that interact by a pairwise two-body 
potential. 

It would be interesting to probe the collective nature of 
fluctuation forces in experiments. So far, only the inter- 
action between two objects, mostly for sphere-plate ge- 
ometries, have been realized experimentally. Among the 
forces studied here, the one between two cylinders and 
a plate, see Fig. [9j might be most feasible in experimen- 
tal studies. To measure the cylinder-plate force, the two 
cylinders need not be separated by vacuum — we expect 
that a similar phenomenon will arise if the cylinders are 
separated by a dielectric spacer layer of fixed thickness. 
This avoids the problem of parallelism one would face 
when measuring the force between two cylinders. Un- 



fortunately, the nonmonotonic effect in Fig. [9] is rather 
small (roughly 0.2%), but it may be possible to increase 
it by further optimization of the geometry and/or the 
material. 

Another important issue is the translational invariance 
of the geometries considered here. The geometries of ex- 
perimental tests will obviously lack this symmetry be- 
yond some length scale. Hence it is important to study 
deviations from the here considered quasi-2D geometries 
due to cylinders and plates of finite size. Then the sur- 
faces have to be treated as compact objects. For this full 
3-dimensional problem, TM and TE modes are longer 
decoupled so that the full electromagnetic vector field 
has to be quantized. This can be done by a recently 
developed multipolc expansion which yields the Casimir 
interaction between objects of arbitrary shape in terms 
of their T-matrices [2"ill22l]. At asymptotically large sep- 
arations between the objects, the interaction is deter- 
mined to leading order by the object's static electric and 
magnetic dipole polarizabilities. For a cylinder of finite 
length L ^> R, the component of the electric polarizabil- 
ity tensor along the cylinder axis scales as L 3 / \og(L/R). 
Hence, the interaction energy between two parallel cylin- 
ders of finite length behaves for separations d L ^ R 

as hc[L 6 /\og 2 (L/R)]d~ 7 . This result shows that for 

large d the interaction is no longer proportional to L or 
the product of the volume of the cylinders as might be 
expected in analogy with the Casimir-Polder interaction 
between spherical particles In fact, for large d the 
interaction amplitude scales as the product of the cubes 
of the largest dimensions of the objects. The cross-over 
between the two extreme cases of infinitely long cylinders 
and asymptotically large separations between finite-sized 
cylinders can be obtained in principle from a multipole 
expansion by including higher order multipolcs. Finally, 
we note that all results given for infinitely long cylinders 
here can be also obtained within the T-matrix approach 
pH . [22| . The latter technique should be particularly use- 
ful to determine whether the same non-monotonic effects 
occur for two spheres next to a metal plate. 
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APPENDIX A: MATRIX ELEMENTS 

Here, we derive M22 and M23 based on Eqs. ([7]) and © 
for TM (Dirichlet) and TE (Neumann) modes in detail. 

We compute M22 for Dirichlet boundary conditions as 
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follows: 



(m\M 22 \m'} = 

G{s 2 L {<t>)- S2 l {4> , );q)e- im <t> +im ''t'' 



2tt 



J Y,e i M-<t>'h j {qr < )K j {qr > y 

3 

(R iq )K m (R iq ). 



-im<f)-\-4m 4>' 



(2n) 2 



(Al) 



The Green's function expansion in terms of modified 
Bessel functions of the first (J) and second kind (K) con- 
tains arguments r< and r> , for which the radius of s^- is 
inserted. Once the angular integrations are carried out 
one obtains in the last line the result of Eq. (fT2]) . 

In order to compute M 22 for Neumann boundary con- 
ditions with the least effort, we make use of the pre- 
vious calculation. According to Eq. ((8]) the directional 
derivatives with respect to the surface normal of cylin- 
ders s?j-(^>) and S2 ((/)') need to be taken inside the inte- 
gral. We simply assume that s 2 (4>) and (</>') are two 
different concentric cylinders with different radii R\ and 
R[, respectively, and i?i < R[ so that the derivatives 
<9 n2 (0)3„ 2 (0') arc taken as d^dR'JjiqR^KjiqR'^. Then, 



the limit R[ — ► R\ is taken to yield Eq. (Q3 

The computation of A/23 is more involved and is carried 
out without a convenient expansion of the Green's func- 
tion. Eqs. (fTU|) and (fTTj) give the surface parametrizations 
which yield 



(m\M 23 \m') = 



2vr 



-jikfs^" (0)— (0') g — imcfy+irn' <j)' 

k 2 + q 2 



d 2 kdcj)d(j)' 
(2tt) 3 



i(k\ R± sin 0+&2 ^1 cos <j>—k\ (R2 sin 0' -\-d) — /C2-R2 cos 0') 



k 2 + q 2 

i*+im'*' d?kd(f>d(f>' 



(A2) 



(2tt) 3 



e -ik lC t / ki +ik2 



k 2 + q 2 y^kj+k 2 



J m (kRi 



k x +ik 2 \ . A 2 k 

X m ,? Jrn'{kR 2 ) — 

V Vfcf + k 2 l 2tt 



where k = \/k\ + k 2 . The limits of integration are, of 
course, to 2tt for the angles <f> and </>' and — oo to oo for 
the k\ and k 2 integrals. Despite the appearance of square 
roots the above expression is analytical in the integration 
variables except for the simple poles due to k\ + k 2 + 
q 2 in the denominator. So, the k\ integration can be 
carried out by contour integration. If d is positive then 
the contour is closed in the lower half plane. The result 



is Eq. (HU): 

(m\M 23 \m') = 
{-i) m+m ' I m {qR 1 )I m ,{qR 2 ) 

(. \ m — m 

V¥+kl+k 2 \ 

I m (qR 1 )I m ,(qR 2 )K m ^ m ,{qd). 



dk2 p-dy/^+kj 



2^fTk, 



(A3) 



M23 for TE (Neumann) modes is easy to obtain. The 
derivatives in Eq. ([8]) can be taken out of the integral 
-they are just derivatives with respect to R\ and R 2 - and 
applied to M 2 3 for TM (Dirichlet) modes. 



APPENDIX B: INVERSE MATRIX AND 
MATRIX ELEMENTS FOR MULTIBODY 
FORCES 

To compute forces according to Eq. ([5]) the matrix M 
or its Fourier transform M need to be inverted. For a 
2x2 block matrix the inverse can be written in terms of 
the inverses of the blocks, 



A B 
C D / 

(A-BD^C)- 1 -A~ 1 B(D - CA~ 1 B)~ 1 
-D- l C(A- BD^C)- 1 (D - CA~ 1 B)~ 1 

(Bl) 



SO, With M np = [Af npj 22, M np ,23; A^np,32, M np ,33] given, 

M^p can be found, then the trace taken, to obtain the 
force. 

To express M np in terms of the matrix elements of 
A/22, A/ 2 ( 2i m, A/23, and M 2 ^ 3R ) we shall refer to their 
parameters in the following way 



M 2 (2,R) = M 2{2 ,R){S,T) 

M 23 = M 23 (d) 

M 2 (3,R) = ^2(3,ft) (d, S, t) ■ 



(B2) 



where the parameters refer to the notation in Eqs. Q13p . 
(|2"3")l . The matrix blocks of Mi p are built up as follows: 



M 1Pi22 = A/22 - s x M 2{2}R) {2H, 



M 1P)23 = M 23 (d) - s x M 2(3tR) (d, 2H, -) 



(B3) 



Clearly, M\ Vj33 is obtained by replacing R\ by i?2 and 
A^ip,32 = ^i P 23- As before, one takes s x = +1 for 
Dirichlet and s x = —1 for Neumann modes. Because 
the reflected cylinders are reflected only once, the minus 
sign is chosen in Eq. (|23[) , as indicated by the minus sign 
in the arguments above. 
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M 2p can be expressed in as follows 

M 2 p,22 = M 22 

„x 



n=Q 

+ M 2 (2,R){- 2H i- 2n (Hi+H 2 ),-) 

OO 

+ Y [^2(2M)(MHl+H 2 ),+) 
n=l 

+ M 2{2 , R) {-2n{H 1 + H 2 ) 1 +) 
M 2 p,23 =M 23 (d) 

OO 

-s x Y, [Mna, R )(d, 2H 1 +2n(Hi+H 3 ), -) 

n=0 

+ M 2{3iR) (d, -2fl 2 -2n(fr 1 +fl 2 ) ) 

OO 

+ Y [M 2{3 , R) (d,2n(H 1 + H 2 ),+) 



2(3, R) 



(d,-2n(H 1 +H 2 ),- 



(B4) 



Again, M2 P ,33 is obtained by replacing R\ by R 2 and 
•^2 P ,32 = M 2p23 . Now cylinders can be reflected an ar- 
bitrary number of times, thus, for an odd number of re- 
flections the minus sign is chosen in Eq. (]23|) and for an 
even number of reflections the plus sign. 



APPENDIX C: ASYMPTOTIC EXPANSION OF 
THE FORCE BETWEEN TWO CYLINDERS IN 
THE PRESENCE OF ONE PLATE 



for N modes are non-diagonal. For H R the non- 
diagonal part can be treated as a small perturbation, 
and the matrix can be inverted pcrturbatively, 



^22 



n=0 



(C4) 



where Af lp j 22oo is the diagonal part of Eqs. (|C2j) . ([C3 
and N is given by 



- ! m{R l\ K m+m ,{2Hq) 



K m ,{Rq)' 



for D modes and 



(C5) 



(C6) 



for N modes. We shall see below that the series of 
Eq. (|C4j) yields a rapidly converging series for the Casimir 
interaction between the cylinders for all sidewall separa- 
tions H > R. Since d is the largest length scale, it is 
useful to set q = u/d in Eq. © since this allows for 
expansions of Bcssel and Hankcl functions that simplify 
further computations. Using Eq. (|C1|) . we obtain for the 
d-dependent part of the energy the series 



E D/N 



hcL 



p OO OO 

/ udu > — Tr 

J° hp 



Awd 2 

Mi P , 32 M^ 22i0O [ Y(-NT ) M lp 



\n=0 



1 P 



2-1 



\n=0 



(C7) 



For simplicity, we consider the case of one sidewall and 
cylinders of equal radius R. The logarithm of the deter- 
minant of the matrix of Eq. (|30[) with n — 1 can be 
expressed as 

logdet(M lp M 1 -| 00 ) 

= -E -Tr(Mr p : 33 Afi P ,32A/r p :22^i P ,23) p (CI) 
p=i P 

+ d-independent terms , 

where we have used log det = Tr log and expanded the 
logarithm. Due to the presence of the sidewall, the "self- 
energy" matrices with elements (see Appendix [Bjl 

(m|Mi Pj 22|m') = (m\Mi p , 33 \m) 

= 6 mm ' I m (Rq)Km(Rq) - I m {Rq)Im' (Rq)K m+m >(2Hq) 

(C2) 

for D modes and 
{m\Mi p , 22 \m') = (m\M 1Pt33 \m') 

= q 2 [Smm'lL{Rq)K' m {Rq) + l' m {Rq)I m , {Rq)K m +m' (2H 

(C3) 



To leading order in 1/d, the integral scales as 1/d 4 . 
At this order it is sufficient to consider the term for 
p = 1 only. Counting the powers of 1/d in the ma- 
trix elements appearing in Eq. (|C7[) shows that at or- 
der 1/d 4 it is sufficient to truncate the matrices at or- 
der m = 1 for D modes and m = 3 for N modes. 
To leading order in 1/d the matrix elements of N are 
~ (R/Hy m \ + \ m I and hence independent of d so that for- 
mally all terms of the series over n have to be included. 
However, since R < H the matrix elements of N n de- 
crease with increasing n, yielding a rapidly converging 
series for the energy. After the expansion in 1 / d the in- 
tegration over u in Eq. (|C7[) can be easily performed for 
N modes since the integrand consists of terms that are 
powers over u and the Hankcl functions K m {u). How- 
ever, for D modes, there are logarithmic corrections of 
the form log(Hu/d) and 1/ \og(Ru/2d). For large d, we 
can make the approximations \og(Hu/d) — » log(H/d) 
and \og{Ru/2d) — » \og(R/2d) since the main contribu- 
tion to the integral comes form u of order unity. With 
these approximations the integration over u can be per- 
formed analytically. By including terms up to n = 5 in 
the series of Eq. (|C7[) . we obtain for the force between 
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the cylinders from N modes (with r = R/H) 



148 i? 4 



55 

1m 1 



625 
9472 ' 



,8 



29 

H 1 

74 

201 

4736 ? 



81 
592' 



.10 



(C8) 



to leading order in 1/d. For D modes, the general form of 
the coefficients of the series in r and its first logarithmic 
correction can be conjectured from the first few terms. 
Rcsummation yields then the closed form expression for 
the force for all r to leading order in 1/d, 



Fd 
L 



he R 4 

77 



G4 



G4 



(4 - 3r 2 ) 5 



+ F 2 (r)log 



( 4 _ r 2)2 

2 {R/2d)- 



\og(R/2d) r 2 (4-r 2 ) 3 
0(log- 3 (R/2d)) +... 



(C9) 



Here i*2(j) = 4/r 4 for r — > so that for H ^> R the force 
is proportional to H 4 /(d 7 \og 2 {R/2d)), which is the same 
as Eq. (|3Tj) up to the different log corrections. Notice 
that the expansion in l/log(i?/2d) of Eq. (|C9j) is formal 
and convergence for all r is not assured in the Dirichlet 



case. The reason for that is related to logarithmic correc- 
tions ~ \og(H/d) = \og{R/2d) + log(2/r) to the matrix 
elements of N so that the expansion coefficients grow 
with H in a way that higher order coefficients involve 
higher powers of \og(H/d). 



The total force is given by the sum Fn + Fd- From this 
large distance expansion we can obtain insight into the 
generation of the non-monotonic behavior of the force. 
The Neumann force expansion is simple to understand, 
Fn decreases in magnitude as the cylinders move away 
from the plate, i.e., r = R/H decreases from 1. Sur- 
prisingly, the expansion of the Dirichlet force without the 
logarithmic terms of Eq. (|C9[) . i.e. at very large cylinder- 
cylinder separations d, indicates that Fd, similar to Fn, 
decreases as the cylinders move away from the plate. But 
ultimately, as H increases, the coefficients of the inverse 
log terms dominate and Fd increases as the cylinders 
move away, albeit at a slower rate compared to Fn, as 
expected from our numerical results and previous con- 
siderations in terms of image charges. If valid, therefore, 
the expansion not only captures the opposing changes in 
Fd and Fn with sidewall separation but also suggests 
that Fd itself has interesting non-monotonic behavior. 
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